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Abstract. 

A renormalization approach has been used in (Eckmann et al 1982) and (Eckmann 
et al 1984) to prove the existence of a universal area-preserving map, a map with 
hyperbolic orbits of all binary periods. The existence of a horseshoe, with positive 
Hausdorff dimension, in its domain was demonstrated in (Gaidashev and Johnson 
2009a). In this paper the coexistence problem is studied, and a computer-aided proof 
is given that no elliptic islands with period less than 20 exist in the domain. It is also 
shown that less than 1.5% of the measure of the domain consists of elliptic islands. 
This is proven by showing that the measure of initial conditions that escape to infinity 
is at least 98.5% of the measure of the domain, and we conjecture that the escaping 
set has full measure. This is highly unexpected, since generically it is believed that for 
conservative systems hyperbolicity and ellipticity coexist. 
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1. Introduction 

The phenomena of coexistence, simuhaneous existence of hyperbohcity, often formulated 
as the existence of positive Lyapunov exponents, and stabihty, in the form of eUiptic 
islands, is believed to be generic in Hamiltonian dynamics (Bunimovich 2008, Gorodetski 
2010, Henon and Heiles 1964, Liverani 2004, Markus and Meyer 1974, Przytycki 
1982, Strelcyn 1991). Very few examples of systems where coexistence has been 
rigorously proved exist; the first such example is given in (Przytycki 1982). At the 
same time very few, if any, systems exist that numerically do not exhibit coexistence, 
as soon as hyperbohcity is present. 

In this paper we demonstrate that the universal area-preserving map associated 
with period doubling does not have any elliptic islands of low periods, and that the 
area of the remaining ones is small. This is a non-trivial result, since for the universal 
map, it is only known that one family of elliptic islands of period 2" have bifurcated 
into hyperbolic periodic orbits. Nothing is known apriori about the periodic orbits with 
period different from 2", and the orbits with period 2" different from those involved 
in the bifurcation sequence. The computations give numerical evidence towards the 
nonexistence of elliptic islands for the universal area-preserving map. Together with our 
resent paper (Gaidashcv and Johnson 2009a), in which we prove that the same map 
has a hyperbolic horseshoe with positive Hausdorff dimension, this gives a candidate for 
an example of a conservative system that does not exhibit coexistence. In fact, these 
results are likely to hold for all infinitely renormalizable maps in a neighborhood of the 
universal one. For area-preserving Henon maps, this contrasts to parameters close to the 
beginning of the period doubling sequence, where the hyperbolic sets are accumulated 
by elliptic islands as proved in (Duarte 2008). 

Our proof is computer-aided and what we actually prove is that there are no elliptic 
islands of period less than 20, and that the area enclosed by the elliptic islands is less 
than 0.049, which is less than 1.5% of the domain. The estimate of the bound on the 
area enclosed by the invariant set scales down like the area of n-ih. generation covers of a 
Cantor set, indicating that the area actually enclosed by elliptic islands is zero. It should 
be noted, however, that the computations are done using the bound on the fixed point 
from (Eckmann et al 1984), and most maps in this neighborhood have elliptic islands. 
Since what we compute are successive approximations of the invariant set, it appears 
that the measure of the invariant set is zero, i.e., that the set of initial conditions that 
escape to infinity has full measure. This is in agreement with the approximations of the 
measure of the invariant set of the conservative Henon family computed in (MacKay 
1982). 

The universal area-preserving map is a map with hyperbolic periodic points of 
all binary periods 2". Its existence, as many other universality phenomena is best 
described in the setting of renormalization. Universality will stand for the independence 
of the quantifiers of the geometry of orbits and bifurcation cascades in families of maps 
of the choice of a particular family. To prove universality one usually introduces a 
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renormalization operator on a functional space, and demonstrates that this operator 
has a hyperbohc fixed point. Universahty was first discovered for the case of period 
doubhng universahty in unimodal maps - the Feigenbaum-CouUet-Tresser universahty 
(Feigenbaum 1978, Feigenbaum 1979, Tresser and CouUet 1978). Proofs appear in e.g. 
(Epstein 1986, Epstein 1989, Lyubich 1999). 

It has been established that the universal behavior in dissipative and 
conservative higher dimensional systems is fundamentally different. The case of the 
dissipative systems is often reducible to the one-dimensional Feigenbaum-CouUet- Tresser 
universahty (Collet et al 1981, de Carvalho et al 2005, Lyubich and Martens 2008, Hazard 
et al 2010). The case of area-preserving maps seems to be very different, and at present 
there is no deep understanding of universality in conservative systems, other than in the 
case of universality for systems near integrability (Koch 2002, Koch 2004, Gaidashev 
2005, Kocic 2005, Khanin et al 2007). 

An infinite period-doubling cascade in families of area-preserving maps was 
observed by several authors in the early 80's (Derrida and Pomeau 1980, Helleman 
1980, Benettin et al 1980, Bountis 1981, Collet et al 1981, MacKay 1993). This 
universality can be explained rigorously if one shows that the renormalization operator 

R[F] = o F o F o A^, (1) 

where Aj;' is some F-dependent coordinate transformation, has a fixed point, and the 
derivative of this operator is hyperbolic at this fixed point. 

It has been argued in (Collet et al 1981) that A^ is a diagonal linear transformation. 
Furthermore, such Ap has been used in (Eckmann et al 1982) and (Eckmann et al 
1984) in a computer assisted proof of existence of a reversible renormalization fixed 
point Fk and hyperbolicity of the operator R. We began the study of the dynamics 
of the renormalization fixed point F* in (Gaidashev and Johnson 2009a, Gaidashev 
and Johnson 2009b, Gaidashev et al 2011). Additional properties of the universality 
phenomena in the area-preserving case arc studied in (Gaidashev and Koch 2011), where 
it is demonstrated that the fixed point F* is very close, in some appropriate sense, to 
an area-preserving Henon-like map 

H*{x, u) = {4>{x) — u,x — 0(0(a:) — u)), 

where solves the following one-dimensional problem of non-Feigenbaum type: 

2 

0(y) = j<f>i<f>i^y)) - y- 

In (Gaidashev and Johnson 2009a) we constructed a hyperbohc set, using the 
method of covering relations (Zghczyhski and Gidea 2004, Zgliczyhski 2009) in rigorous 
computations. To prove that the constructed horseshoe has positive Hausdorff 
dimension wc used the Duarte Distortion Theorem (Duarte 2000, Gorodetski and 
Kaloshin 2009) which enables one to use the distortion of a Cantor set to find bounds 
on the dimension. 

The structure of this paper is as follows: in Section 2 we recall the necessary 
notation and concepts about area-preserving maps, their renormalization theory and the 



No elliptic islands for the universal area-preserving map 



3 



description of the domain of F^,, in Section 3 we summarize the results from (Gaidashev 
and Johnson 2009a) about the hyperbohcity of F*, in Section 4 we state our theorems 
and conjectures, in Section 5 we describe our computer-aided proof, and finally in Section 
6 we describe the results of our computations. 



2. Renormalization of reversible area-preserving maps 

An area-preserving map will mean an exact symplectic diffeomorphism of a subset of 
IR^ onto its image. 

Recall that an area-preserving map can be uniquely specified by its generating 
function S (under the assumption that is a coordinate system): 

^ \ ^ f y 



S,^d,S. (2) 

-Si[x,y)J \S2{x,y); 
Furthermore, we will assume that F is reversible, that is 

ToFoT^F'^, where T{x,u) ^ {x, -u). (3) 

For such maps it follows from (2) that 

Si{y, x) = S2{x, y) = s{x, y), (4) 

and 

^ \ f y 



-s{y,x)J \s{x,y)J' 
It is this little s that will be referred to below as the generating function. We use 
the same normalization as in (Eckmann et al 1984), i.e., s(l, 0) = and d2s{l, 0) = 0.2. 
It follows from (4) that Si is symmetric. If the equation —s{y,x) = u has a unique 
differentiable solution y = y{x, u), then the derivative of such a map F is given by the 
following formula: 



DF{x,u) 



S2{y{x,u),x) 1 

si{y(x,u),x) si(y(x,u),x 
\ S2{y{x,u),x) S2{x,y(x,u) 



(6) 



Si(x,|/(x,M))-S2(x,7/(x,«))fgggg - s,{y{x,u),x) 

The following definition of the renormalization operator acting on generating 
functions is due to the authors of (Eckmann et al 1982) and (Eckmann et al 1984): 



Definition 2.1 



T^EKw[s\{x,y) ^ ji ^s{z{x,y),\y), (7) 

where 

= s{\x, z{x, y)) + s{\y, z{x, y)), 

= s(A, l) + s(0, 1) and ^ = diz{l,0). 
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Definition 2.2 The Banach space of functions s{x,y) — YlTj=o^vi-''^ ~ 0-5)*(y — 0.5)-^, 
analytic on a bi-disk 

\x — 0.5| < p,\y — 0.5| < p, 

for which the norm 

oo 

ikiip ~ \^ij\p 

i,j=0 

is finite, will be referred to as A{p). As{p) will denote its symmetric subspace {s e 
A(p) : si(x,y) = si(y,x)}. 

As we have already mentioned, the following has been proved with the help of a 
computer in (Eckmann et al 1982) and (Eckmann et al 1984): 

Theorem 2.3 There exist a polynomial Sa G As{p) and a ball -B,.(sa) C As{p), 
r = 6.0 X 10"'', p = 1.6, such that the operator TZekw ^■s well-defined, analytic and 
compact on Br- 

Furthermore, its derivative DTZekwIbt has exactly two eigenvalues 6i and 62 of 
modulus larger than 1, while 

spec(DnEKw\Br) \ {Si, C {z eC:\z\ <u < 1}. 

Finally, there is an s* G Br such that 

T^EKw[s*] = s*. 

The scalings A* and /i* corresponding to the fixed point s* satisfy 

A* G [-0.24887681,-0.24887376], 
At* G [0.061107811,0.061112465]. 

Remark 2.4 The radius of the contracting part of the spectrum spec{DTZEKw{s*)) \ 
{61,62} has been estimated in (Eckmann et al 1984) to be u = 0.8. In (Gaidashev et al 
2011) it is proved that v = |A*|. 

It follows from the above theorem that there exist a codimension 2 local stable manifold 

WU^"^) C Br. 

The bound on the fixed point generating function s* will be called s*: 
s* = {s G As{p) : ||s - Sallp < r = 6.0 X 10"^} , 
while the bound on the renormalization fixed point will be referred to as F*: 

Y,. = {F: {x, -s{y, x)) ^ {y, s{x, y)) : s e s*} , 
where Sa is as in Theorem 2.3. 
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We will now summarize the description of the domain of the universal map, F*, from 
(Gaidashev and Johnson 2009a). is defined implicitly by the generating function s*, 
whose domain is given as in (Eckmann et al 1984): 

Vs = {(x, y) e : |x - 0.5| < 1.6, \y - 0.5| < 1.6}. 

To compute the domain of F e F*, we note that its second argument is equal to —s{y, x), 
for some s e s*, see (5). 

Thus, the domain of F, D, is given by: 

V = {{x, u) eC^ : x), s e s*{y, x), \x - 0.5| < 1.6, \y - 0.5| < 1.6}. 

We denote by 

V = {{x,u) eV ■.'^x = '^u = 0}, 

the real slice of V. In (Gaidashev and Johnson 2009a) we constructed a nonempty 
subset of V, which is described in the following lemma and illustrated in Figure 1. 




-1,5 -1 -0,5 0.5 1 1,5 2 2,5 



Figure 1. The real slice of the domain of F* together with its invariant subset. 

Lemma 2.5 There exists a non-empty open set T> (ZM.^, 

VxV g{x eR:\x- 0.5| < 1.6), 

such that for every {x,u) e V and every s e s*{y,x) there exists a unique real solution 
of the equation u — —s{y, x) that satisfies \y — 0.5| < 1.6. 

Clearly, f) (Z T). It also follows that T> contains an open complex neighborhood 
of the set V. The measure of the set V is at least 3.312. With the invariant subset 
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of P, we denote the set of points whose orbits stay in T>. We note that by using the 
renormahzation equation = R[F^], has, for any positive integer k, an analytic 
continuation to the domains 

where we note that each is the rescahng of the subset of a simply connected set 
that is mapped into itself by F*, which is an orientation preserving diffeomorphism. If 
F^,{T>'^~^) n is not connected, we restrict to the connected component containing 
the point po described in Theorem 3.1. This set is simply connected, since D^~^ and 
F^{T>^~^) are. Thus, inductively, for each A;, is simply connected. We put 

oo 

y V^. (9) 

fe=0 

The real slices of the domains are denoted by V^. The three analytic continuations of 
to V^,V^, and V^^ are shown in Figure 2. We note that for all the approximations 
we have computed, F^{V^~^) n V^'^ is connected and V^~^ c V^. Similarly, F~^ has 
analytic continuations to the sets TV^. 




Figure 2. The real slices of the domains and V^^ of three analytic 

continuations of F*. 



3. Review of hyperbolic behavior 

We will now summarize our findings from (Gaidashev and Johnson 2009a) that are 
needed to state the result, the wording has been shghtly modified to better suit the 
present paper: 

Theorem 3.1 The renormalization fixed point has the following properties: 

1) F^ possesses a hyperbolic fixed point po ; 

2) F^ possesses a point which is transversally homoclinic to the fixed point po; 

3) there exists a positive integer n such that for any negative integer k the map F^ 
has a heteroclinic orbit Ok between the periodic points Aj (po) and A^"'"^(po); and for 

nk 

any positive integer k the map F" has a heteroclinic orbit Ok between the periodic 
points Kl{po) and KI~^{pq); 
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Theorem 3.2 = o o admits a hyperbolic set Cq^ C V^, where denotes 
the domain ofG^; 

whose Hausdorff dimension satisfies: 

0.7673 > dimniCG,) > e, 
where e ~ 0.00013 e~^^^^ is strictly positive. 

The transversal homochnic intersection described in Theorem 3.1 is illustrated 
in Figure 3(a), and the two component horseshoe described in Theorem 3.2 is 
illustrated in Figure 3(b). To summarize, in (Gaidashev and Johnson 2009a) we 
prove that the domain, V, of contains: a hyperbolic fixed point located at 
approximately (0.577619, 0), whose stable and unstable manifolds intersect transversally 
at a homoclinic point located at approximately (1.067707,0); a two component 
horseshoe, one component contains the fixed point and the other component contains 
a period three point located at approximately (—0.527155,0); a sequence of period 2" 
points located at approximately A" x (0.577619, 0), with transversal intersections of the 
stable and unstable manifolds of the points with period 2"~^ and 2". 



Figure 3. (a) The homoclinic tangle at the fixed point po. (b) The two components 
horseshoe of G* . 



4. Statement of results 

Theorem 4.1 The map has no elliptic islands of period less than 20 in its domain 
V. 

Most likely, the above result extends to all infinitely renormalizable maps, i.e., maps 
on the renormalization stable manifold, in a neighborhood of F*. This would follow if 
one could prove that such maps are conjugated. Note, however, that the conjugacies 
between infinitely renormalizable maps constructed in (Gaidashev and Johnson 2009b) 
are only defined on the "stable" set, which is the closure of the orbit of (0, 0). It is not 
clear how, and if, these conjugacies can be extended to the entire invariant set. 

From the renormalization equation it follows immediately that if x is a period 2n 
point, then K~^x is a period n point. Hence Theorem 4.1 implies the following corollary. 
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Corollary 4.2 Forn = 1, 2, 3, 4, the analytic continuation of the map to the domain 
has no elliptic islands of period less than [20/2" J in "D". 

Formally, the opposite of Corollary 4.2 is also true: there are no eUiptic islands 
of period 2''n, A; e Z+, n < 20 in A'^V. The invariant subset of A'^V for any A; > 0, 
however, is empty, as can be proved using the methods from Section 5. 

The invariant set, illustrated in Figure 4(a), consists of three regions, from left to 
right: the left component, i.e., the one containing the period three point, of the two 
component horseshoe; the sequence of rescaled periodic points with period 2", and the 
associated multi component horseshoe induced by the transversal intersections of their 
stable and unstable manifolds; the topological horseshoe associated with the transversal 
homoclinic intersection. All of these regions are of a hyperbolic character, and contain, 
according to Theorem 4.1, no elliptic islands of low period. In addition, the invariant 
set also contains the "stable" set from (Gaidashev and Johnson 2009b). 

To give further evidence that there are no elliptic islands, we have computed upper 
bounds on the measure of the invariant subset of V. As is seen from Figure 4(a), most 
of the remaining measure of the cover is around the sequence of rescaled periodic points, 
and their images. We note that the image of the point (0, 0) is close to (1, 0). 

Theorem 4.3 The Lebesgue measure of the invariant subset ofV is less than 1.4594% 
of the measure ofV. 

We will describe in Section 6 that as the measure of the invariant set is estimated 
with higher resolution the estimate scales uniformly with the discretization size, this 
indicates that the measure should go to zero with the discretization size. It also 
shows that, Tip to the discretization size there are no elliptic islands, since once the 
discretization size is smaller than an elliptic island, the estimate of the measure would 
not decrease at the same rate. We therefore make the following conjecture, a similar 
conclusion for the renormalization limit of quadratic area-preserving maps was made 
in (MacKay 1982, MacKay 1993). Note that this statement is saying that the set of 
escaping orbits has full measure in the domain. 

Conjecture 4.4 The Lebesgue measure of the invariant subset ofV is 0. 

Based on Theorems 4.1 and 4.3, and the structure of the invariant set, compared 
to the structure of the periodic orbits, see Figure 4 (a) and (b), we make the following 
conjecture. 

Conjecture 4.5 The analytic continuation of the map to the domain T>°° has no 
elliptic islands in T>°° . 

5. Method 

In this section we describe the algorithms, and their mathematical justifications, that 
are used to prove Theorems 4.1 and 4.3. Our proof uses interval arithmetic, see e.g. 
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-0.4 -0,2 0.2 0.4 0.6 0,8 1 1,2 



Figure 4. (a) The invariant subset of V, and (b) the periodic orbits up to order 19. 
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(Neumaier 1990), to rigorously compute enclosures of the action of the map F*. Based 
on these enclosures we can construct a combinatorial representation of the map using 
graph techniques, i.e., periodic points of the map correspond to cycles in the graph 
of the combinatorial function. This approach is standard for dissipative systems and 
used in e.g. (Gahas 2001, Gahas 2002a, Gahas 2002b, Tucker and Wilczak 2009). 
Since we study a conservative system there are additional difficulties, primarily the 
fact that once a low order cycle is detected, it is impossible to find a trapping region 
for it. Such information is used in (Gahas 2001, Gahas 2002a, Gahas 2002b, Tucker 
and Wilczak 2009) to adaptively reduce the number of nodes in the graph. In the 
conservative case we have to keep all nodes in the graph for all lengths of the cycles, 
there is nothing that prevents the map from having a higher order periodic orbit in a 
previously computed lower order cycle. This significantly increases the complexity of the 
algorithm. Furthermore, the essentially one dimensional discretization method that we 
use, does not lend itself to a nonuniform adaptive partition of the domain, as suggested 
in e.g. (Galias 2002b), instead we use a uniform adaptive discretization of the domain. 
A computational novelty of our paper is the use of generating functions in this setting. 

Our algorithm consists of three steps, which are described in the subsections below. 
First we describe how the reduction to a combinatorial function is performed, second we 
describe how the candidate cycles are contracted, and in the final subsection we describe 
how we check for hyperbolicity. The algorithms described in this section are general, 
and can be used to prove the nonexistence of elliptic islands, up to some period, for any 
arca-prcscrving map. The results of the computation for the universal area-preserving 
map associated with period doubling are given in the next section. 

An implementation (Progs 2010) of the algorithm described below has been made 
using the CAPD library (CAPD 2010) for interval arithmetic operations. 

5.1. Combinatorial representation 

The idea behind the combinatorial representation is simple, partition the domain, 
compute the range of the map restricted to each set in the partition, and check which 
sets in the partition that the image intersects. Once all these ranges are computed, a 
graph is constructed in which each set in the partition is represented by a node, and 
each intersection is represented by an edge. 

Recall our definition of an area-preserving map, it is an exact symplectic reversible 
diffeomorphism of a subset of onto its image. Thus, the action of F^, is described by 
a generating function. Therefore, it is enough to discretizc one of the two dimensions. 
With a physical interpretation of the generating function, we only discretize the position 
variable, not the momentum. 

Consider three consecutive points, {x,Px), {y,Py), {z,Pz), on an orbit of F*, with 
position values x, y, and z. The following holds. 




y 

s{z,y) 




(10) 
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Thus, for any three consecutive points on the trajectory the following equation is 
satisfied: 



Given a pair of points (x, y); (11) is a nonlinear equation in z that can be solved by e.g. 
the interval Newton method. 

To build the graph, we first uniformly discrctizc the projection of the domain on 
the first variable. This is an interval, which we denote by /, the discretization is denoted 
hy U = {Li}f^i. The combinatorial representation is a x A^-graph, where an edge 
between and {Ij,Ik) represents that there exists a trajectory whose projection 

onto the first coordinate visits Jj, Ij, and Ik, in that order. For each pair (Jj, Ij), (11) is 
solved using the interval Newton method. For each set Ik with a non-empty intersection 
with the result, an edge is added to the graph. 

From the graph, all cycles up to some given length are computed by an exhaustive 
depth first search algorithm. All periodic points of the map are contained in a cycle 
whose length divides the period. Due to ovcrestimation and discretization errors, 
however, most cycles are spurious and do not actually contain periodic orbits. All cycles 
detected in the graph are investigated further, as described in the next subsection. 

To decrease the computing time, the search for cycles in the graph is preceded by 
a contraction step, similar to the one in (Galias 2002a, Galias 2002b), where all nodes 
that do not have both inward and outward pointing edges are removed. This procedure 
is performed since an edge that goes to a node which leaves the domain cannot be part 
of an orbit between points in the invariant set. Similarly, a node that is unreachable 
is not in the invariant set. Each time a node is removed, all edges pointing to or from 
it are also removed. This procedure is iterated until no more nodes are removed. The 
remaining nodes represent the invariant set of the domain. In our case, this reduction 
step removes about 90% of the nodes and edges from the graph. 

5.2. Contracting the candidate cycles 

Given a period n, the first step of the algorithm produces a set C„ of sequences of 
n intervals from the partition 



((/|, I2), {I2, 13), (-^n-n -^n)' (-^n> -^i)) graph constructed in the previous section. 

That is, the set C„ is the set of all sequences, P, of intervals, such that there exists a 
cycle in the graph, whose nodes correspond to neighboring pairs of intervals in P, in 
the same order. 

To quickly remove most of the spurious cycles we contract each sequence by solving 
n one dimensional problems, constructed by solving for y in (11), using the interval 
Newton method. Below / = P, for some i. I.e., we iterate: 



Ii = Iin (mid(/i) - (s2(/i-i, li) + S2(/i+i, li)) \s{Ii-i, mid(/i)) + s(/i+i, mid(/i)))) , 



s{x,y) + s{z,y) = 0. 



(11) 




corresponds to a cycle 
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where mid(Jj) denotes the midpoint of the interval Jj. If any intersection in a cycle is 
empty, the entire corresponding cycle is removed. This step removes the vast majority 
of candidate cycles. The remaining cycles are contracted more tightly on the candidate 
periodic orbit using the interval Krawczyk method. We define: 

h) + s(/2, h) 



Z{h,...,Q 



its derivative is given by: 



S{ln-l,ln) + S{ll,ln) 



32 (-Tn 



(12) 



(13) 



,-fl) + S2(-f2,-fl) 



f2(h,l2) + S2{l3, I2) 



Sl{l3, I2) 







si(/„-l,-rn) S2(I„-l,In) + S2{H,In) 



Zeros of (12) correspond to periodic cycles of F*. Note that (13) is a symmetric matrix, 
since Si is a symmetric function. The interval Krawczyk operator is defined as: 

K{I) :- mid(/) - CZ(mid(/)) + {Id - CDZ{I)){I - mid(/)), 

where C is any nonsingular matrix; if possible we use C = (DZ(mid(/)))~^. Similarly 
as with the Newton operator above, we contract the cycles by iterating the following 
intersection: 

7 = 7 n K(I) 

If the intersection of a cycle is empty, it is removed. See (Neumaier 1990) for a thorough 
description of the interval Newton and Krawczyk methods. 

Note that we do not check for existence of periodic orbits, even though the Krawczyk 
method can be used for this. The reason is that we are only interested in proving that 
there are no elliptic cycles, so if we additionally prove that some non-periodic length 
n orbits are hyperbolic that is irrelevant. Furthermore, given the discretization size, a 
candidate cycle might contain several periodic orbits of the same length, and to decrease 
the complexity of the algorithm we want to avoid having to split this interval in several 
subintervals, which would be required if we wanted to prove existence of periodic orbits. 



5.3. Proving hyperbolicity 

Given that the maps that we study are area-preserving, the only possible eigenvalues of 
the derivative at a periodic point are real, or complex and of unit modulus. We compute 
an enclosure of the derivative along each candidate cycle using equation (6), and check 
that the eigenvalues are real and not equal. In the case when we are unable to prove 
whether a cycle is hyperbolic or elliptic, we split the widest enclosure, 7j, of a point in 
the orbit. Note that, 7j C 7j, for some 7j in the cycle, since the original partition is 
contracted. 
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6. Computational results 

In this section we describe the results of the computations using the algorithm presented 
in the previous section on the universal area-preserving map associated with period 
doubhng, F*. These computations constitute the proofs of Theorems 4.1 and 4.3. 
Computations were made on a 32 CPU computer (8 2.2 GHz Athlon Quad Processors) 
with 32Gb of RAM, and on a 3.2 GHz Dual-Core AMD Opteron Processor with 66Gb of 
RAM. The C-|— I- code of the implementation is available for download at (Progs 2010). 

6.1. The invariant set 

The projection of the domain, V, on the first coordinate is / = [—1.1,2.1]. We start 
with a coarse discretization, using 100 subintervals, in order to find oTit whether the 
invariant set is all of / or only a subset. We repeat this procedure a few times until the 
results stabilize. This leaves us with the invariant interval [—0.55,1.13], which is the 
interval that we will study for the remainder of this paper. We set 

7= [-0.55,1.13]. 

To estimate the measure of the invariant set, we compute successive covers of it, as 
described in the previous section. We start with the discretization size 400. Each node, 
(/j, /j), in the graph corresponds to a set, (/j, —s{Ij, li)), in the cover of the invariant set. 
We estimate the area of each box, and sum over all nodes in the graph to get the upper 
estimate of the cover. Given a discretization size, we split each node in two, regenerate 
the edges, and compute the new enclosure of the invariant set. The measure of the 
enclosures at different (nominal) discretization sizes are given in Table 6.1. Each time 
the size of the sides of the covering boxes is decreased by a factor of two, the remaining 
area is decreased by a factor 1.5. Such a uniform scaling is consistent with a cover of a 
Cantor set, giving strong evidence that there are no elliptic islands in the domain. These 
computations also shows that there arc no elliptic islands in the domain that contain 
a ball with diameter larger than 0.00013125, the final discretization size. Since, if such 
islands existed, the improvement would decrease, once this size was passed. Together 
with the bound 3.312 on the measure of V from Section 2, these computations finish 
the proof of Theorem 4.3. 

6.2. Periodic orbits 

Before we start to compute enclosures of the periodic orbits we have to choose the 
(nominal) discretization size of / that will be used for the main computations. After 
some trial and error, it turns out that the most efficient is to choose a discretization which 
is as fine as possible, since this reduces the number of spurious cycles and the number 
of cycles that must be split during the contraction part of the algorithm. Partitioning 
/ into N subintervals means that we will study a (nominally) A'"^-graph, so there is a 
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TiiG/^y/^'I'i rza "I" 1 r\Y\ ( ti r^TTim alt 

Ly±0\jl C UlZjCtulUli I liVjilililCti 1 


Nodes Enclosed Area Improvement 


400 


9366 


0.3742 




800 


25051 


0.2492 


1.50 


1600 


67277 


0.1671 


1.49 


3200 


176304 


0.1097 


1.52 


6400 


459828 


0.07184 


1.53 


12800 


1227236 


0.04834 


1.49 



Table 1. The upper estimates of the measure of the invariant subset of V. The 
nominal discretization size is the number of intervals with uniform size that the interval 
[—0.55,1.13] would be split into if the discretization was non-adaptive, i.e., the first 
coordinate of each box in the cover has length 1.68/discretization. The last column 
shows the improvement of the estimate compared to the previous level. 

memory bound on how large we can take N. The largest (nominal) N that we are able 
to handle conveniently is N — 3200, which is what we use in our proof. I.e. we set 

li = [-0.55 + 0.000525i, -0.55 + 0.000525(i + 1)], for i = 0, 3199. 

This means that, apriori, our graph has 10240000 nodes and up to 32768000000 
edges. Fortunately, our combinatorial function is much smaller. Running the reduction 
procedure, which reduces the graph to the invariant set, at the finest resolution, yields 
a graph with 117883 nodes and 413556 edges. The corresponding invariant set is 
illustrated in Figure 4(a). Ideally, since we study a map, as the discretization size 
goes to zero, the number of edges per node should go to one. With the discretization at 
hand, we get about 3.51 edges per node. 

For each period we locate all the corresponding cycles, contract them, and check 
for hyperbolicity as explained in the previous section. The memory bound part of the 
algorithm is a pre-calculation reduction step that locates which nodes that can be part 
of a cycle of a given length (period 19 requires approximately 100Gb). Wc only search 
for cycles in these reduced graphs. Each reduced graph is much smaller than the graph 
itself, the largest subgraph is the one for period 18, which has 25631 nodes and 69956 
edges. The number of nodes and edges of each subgraph, and the number of cycles, 6„, 
that we found in the graph are given in table 6.2, together with the number of cycles, 
hn, that remained after the contraction procedures. We only checked for hyperbolicity 
on the cycles that remained. The actual number of cycles that we had to compute was 
between 6„ and bn/n, since a cycle of length k appears k times, but we only compute 
it once. This is done by keeping the nodes sorted, and traversing the graph so that no 
nodes with lower index than the starting index is visited. This also implies that we can 
start searching for cycles from all nodes in parallel, which significantly speeds up the 
algorithm. 

The cycles up to order 19 that wc checked for hyperbolicity are illustrated in Figure 
4(b). Note that we do not claim that all of these cycles contain hyperbolic periodic 
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points; what we prove is that all periodic points of order less than 20 are contained 
in these cycles, and that all periodic orbits up to order 19 that are contained in those 
cycles are hyperbolic. This concludes the proof of Theorem 4.1. 



n 


T^nodes 


hedges 

7/ o 


6*1 


L n J 


h„ 


1 


1 


1 


1 


1 


1 


2 


15 


27 


15 


7 


2 


3 


15 


28 


22 


7 


3 


4 


146 


425 


419 


104 


5 


5 


52 


92 


206 


41 


3 


6 


167 


385 


1938 


323 


8 


7 


117 


234 


1807 


258 


3 


8 


1663 


5336 


212547 


26568 


29 


9 


2135 


5811 


265342 


29482 


16 


10 


787 


1838 


188230 


18823 


14 


11 


1060 


2270 


377323 


34302 


21 


12 


5072 


14604 


15861950 


1321829 


187 


13 


1691 


3425 


2123239 


163326 


33 


14 


4275 


9893 


30629187 


2187799 


104 


15 


5106 


11360 


28118492 


1874566 


273 


16 


12823 


35343 


4930683251 


308167703 


80901 


17 


11324 


25850 


253950540 


14938267 


220 


18 


25631 


69956 


24952575180 


1386254176 


13159 


19 


24558 


62806 


3638886917 


191520364 


435 



Table 2. For each period n, the number of nodes and edges refer to the corresponding 
subgraph, 6„ is the number of cycles, including multiplicity, the actual number of cycles 

that arc computed is between fo„ /n and 6„ . /i„ is the number of cycles that remained 
after the contraction and were checked for hyperbolicity. 

The periodic cycles with periods less than 20 arc shown in Figures 5 and 6. To 
illustrate the complexity of the computations, approximate total runtimes for the last 4 
computed periods are given in table 6.2. 



n 


runtime 


16 


3025 


17 


35 


18 


11060 


19 


1495 



Table 3. The total runtimes, in hours, for the proof at different periods. 
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Figure 5. The enclosures of the periodic cycles of of period up to 15. 
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Figure 6. The enclosures of the periodic cycles of F* with periods between 16 and 19. 
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